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Abstract 

O I The possible candidates of 5— wave dibaryons with various strange num- 

bers are studied under the chiral SU{3) quark model. It is shown that there 
are three types of baryon-baryon bound states. The states of the first type 
j^. ■ are called deuteron-like states. If chiral fields can provide enough attraction 

%jj . between interacting baryons, these systems, such as [Sf2 — S*f2] (1^1/2), [SS](o,i), 

[NQ](^2,i/2) would be weakly bound. The states of the second type such as 
[E* A] (0,5/2), [S* A] (3,1/2), [A A] (0,3) and [A A] (3,0) are named as AA-like states. 
Due to the highly symmetric character in orbital space, these systems could 
be relatively deeply bound, but the strong decay modes of composed baryons 
cause the widths of the states much broader. The states of the third type 
are entitled as fifi-like states. Due to the same symmetry character shown in 
the systems of the second type and the only weak decay mode of composed 
baryons, for instance in [r2r2](o,o), or at most one strong decay mode of com- 
posed baryons, for example in [S*r2] (0,1/2), these states are deeply bound states 
with narrow widths. The states of latter two types are most interesting new 
dibaryon states and should be carefully investigated both theoretically and 
experimentally. 



^This work was partly supported by the National Natural Science Foundation of China 
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1. Introduction. 

Dibaryon as a six-quark system lias shown its special place in the investiga- 
tion of medium-energy physics since Jaffe published his prediction of H particle 
in 1977 [0. As is commonly believed, to study a quark system, the Quan- 
tum Chromodynamics (QCD), which governs the strong interaction among 
quarks and gluons, should be employed as an underline theory although its 
non-perturbative behavior is still not quite clear and cannot exactly be solved 
up to now. Jaffe studied the color-magnetic interaction {CMI) of the one- 
gluon-exchange {OGE) potential in the multi-quark system and found that 
CMI shows attractive character in the H particle case. This character com- 
pels six quarks staying in a smaller volume, say less than 0.85fm in radius. 
Thus dibaryon study could provide more information about the short-range 
behavior of QCD, and the existence of dibaryon can directly supply the evi- 
dence of the quark-gluon degrees of freedom in hadrons and hadronic systems. 
However, the reason for forming the baryon-baryon bound state presents great 
complexity. The nonperturbative QCD (NPQCD) effect may seriously affect 
the properties of the dibaryon due to its finite size. The CMI and the in- 
teraction describing the action from physical vacuum should be co-responsible 
for the binding energy of the system. The symmetry structure of the system 
may also play an important role there. In a word, the character of CMI is no 
longer dominant. Sometimes, meson clouds may provide predominant effects. 
Exploring dibaryon may enable us to investigate the short- and medium-range 
NPQCD effects and to find out a practical way to properly treat them. 

In order to reliably study dibaryon, a model that can describe most of short- 
and medium-range NPQCD effects should be employed. In other word, the 
model should have predictive power. It should at least contain two-fold require- 
ments: By using this model the ground state properties of baryons should well 
be fitted, and the experimental data of the nucleon-nucleon (A^ — A^) scat- 
tering and, especially, the empirical data of the nucleon-hyperon (A^ — Y) 
scattering and reaction can reasonably be reproduced in the dynamical calcu- 
lation. When extending this model to dibaryon investigation, no additional 
parameters are required. There are lots of models such as MIT bag model , 
the constituent quark models of various kinds [0, 0, |^, Skyrme model 0, etc. 
Among them, the chiral SU{3) quark model [0, |] is one of the models which 



satisfy requested conditions. In terms of this model, investigating and further 
systematically analyzing possible bound six-quark systems become significant 
and essential. 

Since 70 's, dibaryon has been intensively studied. The most interesting 
dibaryons have been studied are the following: H particle has been theoret- 
ically and experimentally investigated for years. The theoretically predicted 
mass is in a large range [0, H, |, ^, ||, |15], ^, say from 2GeV to 2AGeV. 



The most believed theoretical prediction is around the AA threshold, namely 
around 2.2?)2GeV [^, |^. However, this particle still has not been found in the 
experiment yet. The most recent data showed that the lower limit of the H 



particle mass is about 2.22GeV IJI^. Except the H particle, possible bound 
baryon-baryon systems in the non-strangeness sector were also investigated. 
d* is one of them. There were number of theoretical predictions by using vari- 
ous models, such as the non-relativistic boson-exchange model [|l3, the quark 
cluster model f^, the quark- delo cat ion model |l^, the chiral SU{?>) quark 



model ||20| and etc.. The predicted masses also spread in a wide range. All 



the predictions are below the threshold of the AA channel of 2AQ4GeV, and 
most of them are above the threshold of the strong decay channel, NNirn, 
of 2.154:GeV. d\J^ = 0~,T = 0) is another interesting particle. In the ex- 
periments of double-charge-exchange reactions, it was found that when the 
energy of the incident pion is 50MeV, there exists a resonance with the mass 
of 2.065GeV and the width of 0.5MeV in the processes with a variety of targets 



2^. To explain that phenomenon, one proposed d'. Although there were many 



theoretical attempts |2^ , the theoretical result is still away from the expected 
value. Whether this phenomenon indicates d' is still under discussion. Up to 
now, these three interesting candidates of dibaryons are still not found or con- 
firmed by experiments. It seems that one should go beyond these candidates 
and should search the possible candidates in a wider region, especially the sys- 
tems with multi-strangeness, in terms of a more reliable model such as chiral 
SU{3) quark model. According to this idea, Yu et al. analyzed the six-quark 
system with a simple six-quark cluster configuration ^^. Later, by employing 
the chiral SU{3) quark model, Zhang, Yu et al. studied QQ{S = 0, T = 0) and 
En{S = 1,T = 1/2) PI, g, and Li and Shen explored E*n{S = 0,T = 1/2) 
and EQ — E*Q{S = 1,T = 1/2) IQ. In this paper, we would present a sys- 



tematic study of possible candidates of S— wave baryon-baryon bound states 
in this model. 

The paper is arranged in the following way: The chiral SU{3) quark model 
is briefly introduced in Sect. 2. In Sect. 3, the results calculated by this model 
are given, and the symmetry characters of the system concerned are discussed. 
The effects of chiral-quark field induced interactions on the binding behaviors 
of systems are detailed analyzed in Sect. 4. In Sect. 5, the model parameter 
dependence of the predicted binding energy is further studied. Finally, in 
Sect. 6, the concluding remark is drawn. 

2. Brief introduction of chiral SU{3) quark model. 



Following Georgi's idea [T^, the quark-chiral SU{3) field interaction can 
be written as 

^I = -gchi^L^lpR-lpR^^lpL) , (1) 

with gch being the quark-chiral field coupling constant, ipi and i/jr being the 
quark-left and right spinors, respectively, and 

E = expKA,//], a = 1,2,. ..8 . (2) 

where tt^ is the Goldstone boson field and Xa the Cell Mann matrix of the 
flavor SU(3) group. Generalizing the linear realization of S in the SU{2) case 
to the SU{3) case, one obtains 

8 8 

^ = Y1 "^^Aa - i Yl '^aK , (3) 

a=0 a=0 

and the interaction Lagrangian 

^I = -gch-ip Yl ^^^'^ + ^ Yl ^aKlb W , (4) 

\a=0 a=0 / 

where Aq is a unitary matrix, ctq, .., as the scalar nonet fields and ttq, .., vrg the 
pseudo-scalar nonet fields. Clearly, £/ is invariant under the infinitesimal chiral 
SU{3)l X SU{3)r transformation. Consequently, one can write the interactive 
Hamiltonian as 



Hch = gchF{q^)i) Y ^clK + i Y '^aKl^ -ip ■ 

\a=0 a=0 J 



(5) 



Here we have inserted a form factor F{q^) to describe the chiral field structure 
n]]. As usual, -F(g^) is taken as 



A 



2 \l/2 

CSB 



'""-[tt^-] ' '«> 



and the cut-off mass Acsb indicates the chiral symmetry breaking scale |Tl 



Then, the SU{3) chiral-field-induced quark-quark potentials can be derived 
in the following : 

K.(ni) = -C{gch,m^aAcsB)Xi{m^^,AcsB,rij){Xaii)XaiJ)) 

+ V^Mj), (7) 

2 
TTl 

VTrSrij) = C{gch, m^^,KcsB)jz — - — [^2("^^,, A^sb, rij){ai ■ a^) 

LZTflqiTflqj 

+ (H{m^j,,)-{^fH{kcsBn,)]s.,]{Ui)U3)) , (8) 



with 



2 
-. Tfl 

^LtTlqiTnqj 
2 * 2 

C7(^,,,m,A) = f^- -m , (10) 

47r A"^ — m"^ 

5ij = 3(<Ti-f)((Tj-f)-(a, -aj) , (11) 

and ?7io-a being the mass of the scalar meson and m.,^^ the mass of the pseudo- 
scalar meson. The explicit forms of functions Xi, X2, H and G can be found 
ini- 

As mentioned in |0], the interactions induced by chiral fields describe the 
NPQCD effect in the low-momentum medium- distance range, which is very im- 
portant in explaining the short- and medium-range forces between two baryons. 
To study the baryon structure and baryon-baryon dynamics, one still needs 
an effective one-gluon-exchange interaction V^'^^ which dominates the short- 
range perturbative QCD behavior and a confinement potential ^j™" which 
provides the NPQCD effect in the long distance and confines three quarks to 
a baryon. Then, the total Hamiltonian of the six-quark system can be written 



as 



with 



where 



H = Y: T,-Ta + Y.y^, , (12) 

i i<j 



V^, = Vg^^ + V--f + V^ , (13) 



1 3 1^ 

Vtr^ =-4(A^Ap4-a5(A^Ap , (16) 



and 



Vtt = T.y<^Sn^) + T.y^Sn,) . (17) 

a=0 a=0 

The model parameter should be fixed before calculation. The coupling 
constant Qch is fixed by 



9ch _ /^\2 9nNtv ^^- 

Att ~ ^5' Att Ml 



<? 



:i8) 



and gjf^^/AiT is taken to be the empirical value of about 14. The masses of 
the pseudo-scalar meson m^, m^, m^' and rrix can be chosen as the masses of 
the real vr, 77, rj' and K, and the mass of the scalar meson ctq can be taken as 



mo-o — 625 MeV, according to the relation |T^ 



ml^ = i2mqy + ml . (19) 

In our previous investigation , it was found the N — N and N — Y scatterings 
are not sensitive to the masses of strange chiral fields. In order to reduce the 
numbers of adjustable parameters, rna^{a = 1,..,8) are also taken to be the 
mass of 1]'. The cut-off mass Acsb for various chiral fields is taken to be 

. ^ f 4.2/m-\ for TT, fC and ctq . . 

^^^ \ 5.0fm'\ for r], r^', a', K and e. ^ ^ 



This set of model parameters is called Model I (Set I) which was frequently 
used in our pervious investigations |0, H, |^ . Because the systems studied in this 
paper mostly comprise strange quarks, the strange chiral clouds surrounding 
the baryons become influential. In order to see this effect, we increase the 
masses and corresponding cut-masses of k and e to the values of lAGeV and 
1.5GeV, respectively, which are close to the masses of real mesons with the 



same quantum numbers |25|. This set of parameters is called Model I (Set II). 
When the values of rriu, mg, bu, Qch, "^tt,,! "^o-a ^^^ Acsb are fixed, the one gluon 
exchange coupling constants gu and gs can be determined by mass splittings 
between A and A^, and S and A, respectively, the confinement strengths a^„, 
a^g and a^^ are fixed by the stability conditions of A^, A and S, respectively, 
and the zero point energies a^^, a^^ and a^° are fixed by the masses of A^, A 



and S + n, respectively. 

As is mentioned above, to predict the dibaryon structure, the model should 
be able to reproduce the data of the N-N and Y-N scatterings reasonably. The 
detailed comparison of the theoretical scattering results and the empirical data 
can be found in |^. Here, we only show a typical plot, the cross section of the 
A—p process in Fig.l. In this figure, the solid and dashed curves represent 
the results with Sets I and II, respectively. It is shown that both curves are 
consistent with the experimental data. After confirming the model, we use the 
same set of parameters to study the dibaryon system. 

All the model parameters in Sets I and II are tabulated in Table 1, respec- 
tively. 

Table 1. Model parameters 





Set I 


Set II 


mu{MeV) 


313 


313 


ms{MeV) 


470 


470 


buifm) 


0.505 


0.505 


9u 


0.936 


0.936 


9s 


0.924 


0.781 


al^{Mev/fm^)/at{Mev) 


54.3/-47.7 


55.0/-48.9 


aUMev/fm')/af,{Mev) 


65.8/-41.7 


66.5/-50.6 


aUMev/fm^)/at{Mev) 


103.0/-50.6 


115.4/ -73.7 


m^{fm-^)/K{fm-^) 


0.7/4.2 


0.7/4.2 


mK{fm-')/AK{fm-') 


2.51/4.2 


2.51/4.2 


m^(/m-i)/A^(/m-i) 


2.78/5.0 


2.78/5.0 


m^/(/m-i)/A^/(/m-i) 


4.85/5.0 


4.85/5.0 


m^,{fm-')/A^,,{fm-') 


3.17/4.2 


3.17/4.2 


m^/(/m-i)/A^,(/m"i) 


4.85/5.0 


4.85/5.0 


m.(/m-i)/A«(/m-i) 


4.85/5.0 


7.09/7.6 


m,(/m-i)/A,(/m-i) 


4.85/5.0 


7.09/7.6 



The binding energy of the baryon-baryon system is dynamically solved by 
using the Resonating Group Method (RGM). In this method, the trial wave 
function of the six-quark system can be written as 



^ = ^[0A(6,6)0B(6,C5)Xre/(^)XCM(^CM)]sr, 



(21) 



where 4>a{b) denotes the antisymmetrized wave function of the baryon cluster 
A(B), Xrei{R) the trial wave function of the relative motion between interacting 
clusters A and B, Xcm{Rcm) the wave function of the motion of the total 
center of mass, and C,i the Jacobi coordinate with i = 1 and 2 for cluster A and 
i = 4 and 5 for cluster B, respectively. The symbol A describes the operation 
of the antisymmetrization between quarks in two interacting clusters. This 
operator can be read as 

A = NY,epP, 
p 

where P is an operator which permutes quarks of cluster A and quarks of 

cluster B, ep = 1(— 1) when P is an even(odd) permutation and N is the 

normalization factor. Considering the permutation symmetry, A can also be 

written as 



A 



N'{l-J2 P^^"^ 



V 



f22) 



jeB 



where P°f denotes the permutation operation carried out between the i — th 
and j — th quarks in the orbit, spin, flavor and color spaces, simultaneously, 
and again N' represents the normalization constant. Calculating the expec- 
tation value of the Hamiltonian operator on the trial wave function in which 
the unknown Xrei{R) is expanded in terms of well-known bases, one deduces a 
secular equation 

J2{H^J - EN,,)c, =0, (z = 1, 2, ■ ■ ■ , n), (23) 

with Hij and Nij being the Hamiltonian and normalization matrix elements, 
respectively, E the eigenvalue and Cj the corresponding eigenfunction, namely 
the expansion coefficients of Xrei{R)- Solving this equation for Cj, one ob- 
tains the binding energy and the corresponding wave function of the six-quark 
system. The detailed method can be found in our previous paper [^^ Q. 



3. Symmetry character discussion. 

As a comprehensive survey, there are two crucial physical factors which 
resolve whether a two-baryon system is bound. One is the symmetry property 
of the system, namely the characteristics of quark exchanges between baryons, 
and the other is the interaction including both the direct and exchange com- 
ponents between quarks, and eventually between baryons. In this section, we 
analyze the symmetry property of the system according to the character of the 
matrix element P^i'^, where the superscript sfc denotes the P operator acts 
within the spin — flavor — color space only and the subscript 36 represents 
the exchange operation is between the 3-rd and 6-th quarks. Then, in section 

4, we discuss the effects of various chiral-quark field interactions by employing 
several models. 

As is pointed out in Ref.[21], the matrix element (^'^■^^) is an important 
measure of the action of the Pauli principle in the two-baryon state. This mea- 
sure specifies the symmetry character of the state of the system. According to 
the symmetry characters of systems, namely the mentioned matrix elements, 
the two-baryon systems concerned can generally be divided into three classes. 
In the first class, QiP^i'^) ~ 1, namely, (A^^^) ~ 0. The Pauli blocking effect 
between interacting baryons are incrediblely serious so that the two-baryon 
5— wave state with [6]^, the [6] symmetry in the orbital space, is almost a for- 



bidden state. Namely, it is very hard to form a bound state. The state in the 
second class has the property of 9(P3g'^) ~ 0, or (A^-^^) ~ 1. In this class, the 
Pauli blocking effect between interacting baryons are very small so that the 
exchange effect between quarks which are located in different baryons, respec- 
tively, becomes negligible and these two baryons are relatively independent 
with each other. As a result, the meson-exchange effect may play a dominant 
role in binding. If the inter-baryon interaction shows an attractive feature with 
a large enough strength, the system would be bound. This kind of system may 
also be deduced in terms of a model in the baryon-meson degrees of freedom. 
The state in the third class possesses a feature of 9(P3g ^) ~ —1, or (A^-^'^) ~ 2. 
The inter-baryon quark-exchange feature of this kind would be enormously 



beneficial to form a state with the [6],. symmetry |21]. If the inter-baryon 
interaction demonstrates the attractive character with certain strength, it is 
possible to form not only a bound state, but also a dibaryon with a relative 
smaller size in radius. We present the resultant binding energies and the cor- 
responding root-mean-square radii {RMS) of the 5— states which have various 
strange number S and belong to the second and third classes in the following 
subsections. 

3.1. The systems in which the expectation values of P^^ operator 
are close to zero. 

In the deuteron case, {P^q^) = —1/81, it is a typical case of the second 
class. We collect some systems which have the same symmetry characteristics 
as deuteron in this subsection. The binding energies, Et, and corresponding 
RMS, TZ, of these systems are tabulated in Table 2. 

Table 2. Binding energy, E),, and corresponding RMS, TZ, for the systems with (^35) ~ 0. The 
units for Ef, and RMS are in MeV and fm, respectively. 



s 


System 


(^3^) 


Model I (Set I) 


Model I (set II) 

Eb l/n 





NN(S = 1,T = Q) 
deutoron 


-1/81 


0.2 //1. 63 


-0.7 //I 


68 


-1 


NA 
(S = 0,T = l/2) 


NA{S = 0,T = 1/2) 





-6.6 //1. 63 


-7.7 //I 


69 


AfA- ALECS' = 0,T = 1/2) 




-5.7 //1. 59 


-7.2 //I 


66 


NA 
(5 = l,T = l/2) 


AfA(S' = l.T = 1/2) 





-7.7 //1. 68 


-7.7 //I 


68 


A'A-ArE(5 = l,T = l/2) 




-6.6 //1.63 


-6.7 1/1 


64 


Ni:(S = 0,T = 3/2) 


-1/81 


-5.1 //1. 58 


-5.3 //I 


59 


-2 


AA 
(5 = 0,T = 0) 


AA(5 = 0,T = 0) 





-4.8 //1. 54 


-5.6 //I 


58 


AA-A^H-SS(5 = 0,T = 0) 
H particle 




-2.0 //1. 15 


8.2 //0.91 


-3 


AH(5 = 1,T = 1/2) 


4/81 


-8.1 //1. 74 


-7.7 //1.72 


7VQ(5 = 2,T = l/2) 





3.5 //1.18 


12.7 //0.98 


AQ{S = 3,T = 3/2) 





4.4 //1. 15 


14.2 //0.96 


-4 


H5(5 = 0,T = 1) 


-1/81 


4.1 //1. 17 


0.4 //1. 30 


-5 


(5 = l,T = l/2) 


Hn(5 = l,T = l/2) 


1/81 


9.5 //1. 02 


4.2 //1. 14 


sn-s*n{s = i,T = i/2) 




32.9 //0.78 


32.6 //0.77 



The data in Table 2 shows that the deuteron is weakly bound in the Model I 
(Set I) case, which indicates that the chiral SU{?)) quark model, in principle, 
can reasonably describe the structure of deuteron^. It is also seen that in the 
single A^A channel and the coupled A^A — NH channel with S = O(orl) and 
T = 1/2 cases and the single iVS channel with S = and T = 3/2 case, no 
matter which set of model parameters is employed, the systems are not bound. 
These are in agreement with experiments. The H particle is also not bound in 
Set I but weakly bound in Set II. The resultant mass of H is close to the AA 
threshold in both Set I and Set II, and this feature is consistent with the recent 



finding in experiments |T^. These results further convince ourselves that the 
chiral SU{3) quark model is reasonable and reliable in the bound-state study. 

It should be noted that (Ps^g ) of these systems appoximately being zero 
means that the symmetry structure of this kind makes the quark-exchange 
effect less important, and consequently, the contribution from the kinetic en- 
ergy term shows relatively repulsive nature to the kinetic energy of the relative 
motion between two well-separated interactive baryons (see Appendix), which 
makes interacting baryons apart. Therefore, very similar to deuteron, whether 
the two-baryon system is bound depends on the feature of the interaction be- 
tween interacting baryons, especially that caused by chiral fields, namely the 
overall characteristics of the short- and medium-range NPQCD effects, to a 

^ Of course, it is easy to fine-tune the mass of a^ so that the best agreement between 
the calculated and experimental binding energy of deuteron can be achieved. However, it is 
not necessary to give a very accurate mass of predicted dibaryon, we would relinquish this 
adjustment. 
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considerable extent. If the characteristics is attractive in nature, the system 
would be bound like deuteron, and we call it as a deuteron-like system. If 
it shows weak attraction or even repulsive feature, the system would not be 
bound anymore. It is also noticed that the OGE interaction provides repul- 
sion in the deuteron, iVA, A^S, SS(5 = 0,T = 1) and ~.VL{S = 1,T = 1/2) 
systems. To form a bound state, a strong enough attractive interaction by 
the meson exchange must be requested. For instance, the vr exchange, espe- 
cially the tensor force, causes the weakly binding of deuteron and the relative 
strong attraction from the coupling between the chiral fields and quarks makes 
SS(S' = 0,T = 1) and 30,(3 = 1,T = 1/2) bound. Then in the H particle 
case, although OGE interaction provides attraction, due to the repulsive ef- 
fect from the kinetic energy term, the only attraction from OGE potential 
is not strong enough to cause binding. The attractive feature of the chiral 
field would be helpful to form a weakly bound H. But it is model parame- 
ter dependent. In the NQ{S = 2,T = 1/2) and AQ{S = 3,T = 3/2) cases, 
OGE contributes nothing, the weakly bound behavior of these systems fully 
depends on the attractive features of the chiral fields. Moreover, it is notewor- 
thy that very similar to the deuteron, the AA — NE — SS(S' = 0, T = 0) and 
En - E*n{S = 1,T = 1/2) systems are bounder than the AA{S = 0,T = 0) 
and EQ{S = 1,T = 1/2) systems due to the channel coupling effect. 

In summary, we predict that Nn{S = 2,T = 1/2), An{S = 3,T = 3/2), 
EE{S = 0,T = 1) and EQ{S = 1,T = 1/2) are weakly bound baryon-baryon 
states. The mass of the H particle is around the AA threshold. EQ — E*Q{S = 
1,T = 1/2) system is a bound state with a relative large binding energy due 
to the strong channel coupling. 

3.2. The systems in which the expectation values of P^q operator 
are close to —1/9. 

The systems in which (-Pjg ) ~ —1/9 are collected in this category. The 
binding energies, Eb, and corresponding root-mean-square radii, TZ, of these 
systems are tabulated in Table 3. 

Table 3. Binding energy, E),, and corresponding RMS, TZ, for the systems with (^35) = —1/9. 
The units for Ej, and RMS are in MeV and fm, respectively. 
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5 


System 


(^3^) 


Model I (Set I) 


Model I (set II) 





AA{S' = 3,T = 0) (d*) 


-1/9 


22.2 //1. 01 


18.5 //1. 05 


AA(5 = 0,T = 3) 


-1/9 


16.0 //1. 10 


13.5 //1. 14 


-1 


i;*A(S = 0,T = 5/2) 


-1/9 


24.6 //0.99 


19.0 //1. 04 


S*A(S' = 3,T= 1/2) 


-1/9 


25.9 //0.95 


29.3 //0.93 


-5 


~*n{S = 0,T = 1/2) 


-1/9 


92.4 //0.71 


76.5 //0.72 


-6 


f7n(5 = 0,T = 0) 


-1/9 


116.1 //0.66 


98.5 //0.67 



/^From this table, one sees that all the systems in this category have the fea- 
ture of (-P36 ) = — 1/9. It indicates that the system has relatively higher 
anti-symmetry in the spin- flavor-color space. As a consequence, the contri- 
bution of the kinetic energy to the binding energy plays a relatively attrac- 
tive role in comparision with the kinetic energy in two independ baryons |^ 



(see Appendix). This characteristics would bring six quarks closer. If the 
chiral fields can additionally provide attraction between interacting baryons, 
the deeply bound state may be established, such as f2f2(S' = 0,T = 0) and 
'E*Q{S = 0,T = 1/2). Even the contributions of chiral fields are mild, the 
system with the symmetry structure of this kind may still have the binding 
energy of several tens MeV. Therefore, in any case, the system in this cat- 
egory would be a bound state. However, it should be noticed that in most 
states here, such as AA and S*A, both composed baryons have strong decay 
modes. Therefore, if the AA and E*A are not bound deeply enough, namely 
their masses are not smaller than the thresholds of NNtctt and NAtctt, re- 
spectively, these four states should be bound states with broad widths. Only 
^2^2(5 = 0,T = 0) and 'E*Q{S = 0,T = 1/2) are the most interesting systems. 
Both f2's in ^2^2(5* = 0, T = 0) can decay only through the weak mode, so that 
QQ{S = 0,T = 0) is a bound state with a narrow width. It is also possible 
that the mass of 2*^2(5 = 0,T = 1/2) is smaller than the threshold of SQir 
(with Set I), then this state could also be a narrow width bound state. 

4. Effects of interactions induced by chiral-quark field couplings. 

In this section, we would demonstrate another factor which dominates the 
binding behavior of the system concerned, namely the interactions caused by 
the chiral-quark field couplings. 

The characters of chiral fiels are different case by case, and the importance 
and sensitivity of these fields in the bound-state problem of six-quark systems 
also dissimilar. In order to analyze the effect of the chiral field on the binding 
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energy of the baryon-baryon system, we build up other two models. In Model 
II, the 77, K, 1], Vj and ctq are considered. This is so called extended chiral 
SU{2) quark model. The Model III is called chiral SU(2) quark model, in 
which only vr and ctq fields are remained. 

We arrange all concerned systems which are possibly bound into following 
types: 

4.1. Deuteron-Iike systems. 

To convince readers the the importance of the interaction for binding, we 
present the contributions of various terms, such as kinetic energy, OGE, pseu- 
doscalar meson exchange and scalar mesons exchange terms to the total bind- 
ing energy for the typical case SS(o,i) in table 4. 

Table 4. Contributions of various terms to binding energy for HH/q i), the unit for energy is 
in MeV. 





kine. 


OGE 


TT 


K 


V 


V 


a 


a' 


K 


e 


Model I Set 1 


-17.3 


-15.2 


-1.2 


-1.7 


-1.4 


-0.8 


37.7 


1.9 


-1.2 


6.2 


Model II 


-12.0 


-5.5 


-0.8 


-1.1 


-0.8 


-0.5 


26.0 


- 


- 


- 


Model III 


-16.1 


-15.8 


-1.4 


- 


- 


- 


40.5 


- 


- 


- 



From the table, we sees that although the kinetic energy, OGE and pseu- 
doscalar mesons provide the repulsive feature, the strong attractive effect from 
scalar mesons, especially a meson, would compensate the repulsion and make 
the SS(o,i) system weakly bound. 

We then demonstrate the binding energies and the corresponding root- 
mean-square-radii of the bound systems in Sect. 3.1 with Models I (set I), II 
and III, respectively, in Table 5, except the H particle which was carefully 
studied in our previous paper |p. 

Table 5. Binding energy, Ei,, and corresponding RMS, TZ, for deuteron-like systems in various 
chiral quark models. The units for E), and RMS are in MeV and fm, respectively. The system is denoted 
by symbol [BiB2]is,T) with Bi, B2, S and T being baryons 1 and 2 and the total spin and isospin of the 
system, respectively. 



Model 


5 = -5 


5 = -4 


5 = -3 


5 = 


pf7-S-n](i,i/2) 

Et //n 


|2SJ(o,l) 

Et //n 


lA^^I (2,1/2) 

Eb //7^ 


iAH] (3,3/2) 

Et //n 


Et lin 


I (Set I) 


32.9 //0.78 


4.1 //1. 17 


3.5 //1. 18 


4.4 //1. 15 


0.2 //1. 63 


II 


29.3 //0.79 


-0.5 //1. 33 


31.8 //0.81 


34.3 //0.80 


2.1 //1. 52 


III 


17.6 //0.86 


3.1 //1. 18 


49.5 //0.74 


49.6 //0.74 


4.4 //1. 41 



As mentioned in the preceding section, the binding energies of the systems 
in this type are sensitive to the contributions offered by the chiral field. This 
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character is very similar to that of deuteron and can exphcitly be seen in this 
table. In most cases in this type, say the systems with low strange number, 
when one switches the model from Model II to Model I (Set I) or even from 
Model III to model II, namely the strange clouds are taken into account step 
by step, additional strange-cloud-caused interactions would make the overall 
inter-baryon interaction less attractive in nature, and consequently the systems 
would become less bounder. In a word, these systems are bounder in the 
chiral SU{2) quark model than in the chiral SU{3) quark model. But, in high 
strangeness systems, additionally taking strange clouds into account would 
benefit the binding. 

The channel sketches and the predicted mass ranges of these states with 
the different combinations of chiral fields, which are denoted by shaded areas, 
are plotted in Figs.2(a)-(d) for the [EQ — S*i7] (1^1/2), [2S](o,i), [A^fi] (2,1/2) and 
[Afi] (3^3/2) systems, respectively. 

4.2. AA- and r2r2-like systems 

Now, we analyze the role or impoatance of the interaction in the system 
in the third class or in subsection 3.2.. As is mentioned above, except the 
symmtry character (^3*6 ) ~ — 1/9 of the system, the effect of the interaction is 
dominantly responsible for binding. A detailed study of the binding energies of 
the systems in this category show that, different with the contribution provided 
by the interaction in the second class, the contribution from exchange term of 
the interaction is substantially large and gernerally has two different character 
for various states in the class. Considering this difference together with the 
different decay mode, we further distinguish the states in this class into AA- 
like states and f2f2-like states. In AA-like states, where only weak decays 
exist, the direct and exchange terms of interactions play comparable roles in 
forming dibaryon and the binding energies of these states are only a few tens 
of MeV. As a specific example, the calculated result showes that in AA(o,3), 
the exchange term of the interaction induced by a meson provides a binding 
energy of 17.0 MeV, which is comparable to the direct contribution of 16.0 
MeV. While in f2f2-like states, where at most one strong decay mode exists, 
the exchange terms of interactions play much more important roles than the 
direct terms do and the binding energies of these states can reach nearly one 
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hundred MeV. For instance, in f2f2(o,o) , the contribution of the exchange term 
of the a interaction is 80.6MeV, which is nearly two times larger than that of 
the direct contribution of 44.3 MeV. We would investigate these two types of 
states in detail in the following. 

4. 2. a. AA-like systems. 

We firstly single out AA(o,3) as a representative to decompose the binding 
energy of the state in this class. The results are tabulated in Table 6. 

Table 6. Contributions of various terms to binding energy for AA^ogj, the unit for energy is 
in MeV. 





kinc. 


OGE 


7r 


K 


V 


n' 


O" 


cr' 


K 


€ 


Model I Set 1 


8.0 


-26.7 


-15.0 


0.0 


-2.3 


-0.4 


33.0 


14.5 


0.0 


4.8 


Model II 


6.2 


-12.7 


-10.0 


0.0 


-1.2 


-0.2 


24.1 


- 


- 


- 


Model III 


12.5 


-17.8 


-17.8 


- 


- 


- 


36.4 


- 


- 


- 



It is seen that only OGE and pseudosclar mesons provide repulsive contri- 
bution, while both kinetic energy and the scalar mesons devote their attractive 
contribution to binding. Moreover, it is shown that again the potential terms, 
especially the terms provided by the scalar meson, dominate the binding be- 
havior. 

The binding energies and the corresponding root-mean-square-radii for this 
type of states with Models I (set I), II and III are tabulated in Table 7, re- 
spectively. 



Table 7. Binding energy, Ei,, and corresponding RMS, TZ, for AA-like systems in various chiral 
quark models. The units for Ef, and RMS are in MeV and /m, respectively. The system is denoted by 
symbol [B1B2] (s,t) with Bi, B2, S and T being baryons 1 and 2 and the total spin and isospin of the 
system, respectively. 



Model 


5 = -1 


5 = 1 


IS* AJ (0.5/2) 

Et //n 


lS'A](3,i/2) 

Eb //n 


iAA](o,3) 

Et //n 


Et //n 


I (Set I) 


24.6 //0.99 


25.9 //0.95 


16.0 //1. 10 


22.2 //1. 01 


II 


11.7 //1. 13 


78.0 //0.79 


6.3 //1. 25 


64.8 //0.84 


III 


21.5 //0.99 


76.3 //0.80 


13.2 //1. 11 


62.7 //0.86 



It is shown that the binding characters vary with the spin structures of 
systems. When one turns the model from Model I (Set I) to Model II, the 
binding energies of the S* = systems become smaller and those of the S = 3 
systems become larger. However, if the model is switched from Model II to 
Model III, in comparison with the preceding case the binding energies for 
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different spin systems change in opposite directions. Nevertheless, substituting 
Model I (Set I) by Model III would cause the binding energy decreasing in the 
S* = system and increasing in the S* = 3 system. Precisely, the systems with 
S = and S* = 3 are more bound in the chiral SU{3) quark model and the 
extended chiral SU{2) quark model, respectively. Moreover, because CMI 
shows repulsive feature in the S* = system and attractive nature in the S* = 3 
system, the binding energies of the S* = 3 systems are generally higher than 
those of the S* = systems. 

We also draw channel sketches and predicted mass ranges of [S* A] (0,5/2), 
[S* A] (3 1/2), [A A] (0,3), [A A] (3,0) states with different combination of chiral 
fields, which are signified by shaded areas, in Fig. 2(e) and Fig. 2(f), respec- 
tively. Since in the S*A cases, the predicted masses are above the threshold 
of the strong decay channel ANtctt, and in the AA cases, the resultant masses 
are above the threshold of the strong decay channel NNtttt, these bound states 
would have large widths. It might not be easy to detect in experiments. 

It is worth to mention that because the Hamiltonian for the octet and 
decuplet baryons are in the same form and the calculated masses of these 
baryons can well-fit the data, the RGM calculations for the systems involving 
A are valid. Of course, whether this model treatment agrees with reality should 
be verified by the experiment in future. 

4.2.b. i7r2-Iike systems. 

As what we did for deuteron- and AA-like states, we choose ^2^2(0,0) as a 
specific example to show the contributions of various Hamiltonian terms. 

Table 8. Contributions of various terms to binding energy for f2f2(oo); the unit for energy is 
in MeV. 





kinc. 


OGE 


n 


K 


V 


V 


a 


a' 


K 


e 


Model I Set 1 
Model II 
Model III 


8.6 
12.3 

12.7 


-82.7 
-33.3 

-58.5 


0.0 
0.0 
0.0 


0.0 
0.0 


-7.6 
-6.5 


-15.0 
-12.5 


124.9 

112.5 

97.7 


0.0 


0.0 


87.7 



The large binding energy of ^2^2(0,0) in model I and the binding behaviors of 
the system in model II and III can clearly be undestood by the detailed con- 
tributions of interactions in table 8. Apparently, the kinetic energy contribute 
attractively to the binding energy and the a meson (in Model I , also e meson) 
provides strong attractive contribution and dominates the binding behavior. 
OGE, f] and rj' are repulsive behavior contributors. 
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We tabulate the binding energies and the corresponding root-mean-square- 
radii for this type of systems with Models I (set I), II and III in Table 9, 
respectively. 

Table 9. Binding energy, Ei,, and corresponding RMS, TZ, for HQ-like systems in various chiral 
quark models. The units for Ei, and RMS arc in MeV and fm, respectively. The system is denoted by 
symbol [i?ii?2](s,T) with Bi, B2, S and T being baryons 1 and 2 and the total spin and isospin of the 
system, respectively. 



Model 


5 = -6 


5 = -5 


If^f^l(O.O) 

Ei l/n 


[s*r2](o_i/2) 
Ei, //n 


I (Set I) 


116.1 //0.66 


92.4 //0.71 


II 


74.2 //0.69 


58.8 //0.75 


III 


53.5 //0.73 


52.1 //0.76 



It is seen that no matter in which model, these two states are bound states, 
or dibaryons. The CMI of the system in this type shows repulsive feature 
which offers disadvantage to binding. However, the benefit from the symmetry 
structure discussed in Sect. 3. 2 somehow overcomes this disadvantage. On the 
other hand, because these systems are high-strangeness systems, the attractive 
characteristics from the a^- field-caused interaction and strange-cloud-induced 
interactions make the systems deeply bound and corresponding TZs less than 
0.76 fm. In other word, by using the chiral SU{3) quark model, these systems 
could be even bounder. More detailed analysis for [r2r2](o,o) is given in Ref. [21]. 

The channel sketches and the predicted masses of the [r2r2](o,o) and [S*f2] (0,1/2) 
states with different combinations of chiral fields, which are presented through 
shaded areas, are shown in Fig. 2(g) and Fig.2(h), respectively. Since Q does 
not have strong decay mode, the predicted [r2r2](o,o) bound state should have 
narrow width. However, in the [S*r2] (0,1/2) system, because S* has the strong 
decay mode and the predicted masses of S*r2 by using different models stride 
over the threshold of the strong decay channel SQtt, whether the width of the 
bound state is narrow depends on the model employed. With the chiral SU{3) 
model, which, we believe, is a right model for the high strangeness system, the 
width should be narrow too. 

Furthermore, we would emphasize that the purpose of this section is to 
trace the influences of various chiral flelds on the calculated results through 
employing different models such as Models II and HI. In our consideration, the 
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result with Model I is the most reliable one, because investigating a system 
with strange quarks without considering strange clouds is unreasonable. 

5. Model parameter dependence 

To ensure the reliability of our predictions, various kinds of model parameter- 
dependencies are also investigated. 

It is clear that the variation of 6„ would change the size of the baryon as well 
as the values of model parameters via the conditions mentioned in Sect. 2, and 
the variation of rris would change the strength of the OGE potential, Qs- We 
plot Eh with respect to 6„ with different rris for six most interesting systems, 
where (-Pse) = —1/9) ^^ Fig-3. From this figure, one finds that the binding 
energy decreases with increasing values of 6„. When 6„ is altered to a large 
value, the higher spin the system has, the milder the binding energy changes, 
while with increasing the mass of strange quark, the more strange quarks the 
system contains, the larger binding energy the system has. As long as the 
value stays in a reasonable region, say from 0.5/m to 0.55/m, these systems 
remain the bound nature. 

The mo--dependence is also studied, because the a induced interaction is 
always a dominant piece that responses for the binding behavior for the two- 
baryon system. The results with different values of m„ are tabulated in Table 

7. 



Table 10. The mcr —dependence of the binding energy Ei,. The units of Ei, and m^ are in MeV. 
The system is denoted by symbol [BiB2](s,T) with Bi, B2, S and T being baryons 1 and 2 and the total 
spin and isospin of the system, respectively. 



m^ 


550 


600 


625 


650 


700 


Eb 


^(0,0) 


140.6 


121.9 


116.1 


107.1 


94.3 


=*f^(0,l/2) 


116.2 


99.5 


92.4 


86.5 


75.4 


=^(1,1/2) 


21.3 


11.2 


9.5 


4.7 


0.5 


JV^(2,l/2) 


17.7 


7.1 


3.5 


0.9 


-2.8 


2*^(0,5/2) 


36.5 


28.1 


24.6 


21.9 


17.0 


^ ^(3,1/2) 


40.8 


30.2 


25.9 


22.5 


16.4 


AA(o,3) 


25.5 


18.7 


16.0 


13.9 


10.2 


AA(3.o) 


35.0 


25.8 


22.2 


19.3 


14.2 


7ViV(i,o) {d) 


5.1 


1.6 


0.2 


-0.5 


-2.1 



lYiom. this table, one sees that as long as rrio- is in a reasonable region, 
namely the existent experimental data can unifiably and reasonably be ex- 
plained, most systems studied would not change their binding behaviors, ex- 



cept those whose binding energies are marginal, for instance, when rricr = 
700MeV, NQ becomes unbound. 

The cutoff mass dependence is also investigated. Our result shows that 
if one increases all the cutoff masses to 1.2GeV, the binding energies of the 
systems concerned would slightly arisen. The binding behaviors are relative 
stable. 

In a word, as long as the model parameters are selected in reasonable 
regions and all the experimental data, such as A^A^ and NY scattering data, 
the ground state masses of baryons, binding features of well-known baryon- 
baryon states, and etc., can be basically reproduced, the concerned six two- 
baryon systems with {P^q'^) = — 1/9 should be bound, especially, f2f2(o,o) and 
S*r2(o,i/2) should be deeply bound dibaryons. 



6. Concluding remarks. 

The possible candidates of S*— wave baryon-baryon bound states with var- 
ious strange numbers are systematically studied by using the chiral quark 
model. In terms of the chiral SU{3) quark model, which are believed to be one 
of the best models in the quark-gluon degrees of freedom in well reproducing as 
much experimental data as possible, we predict the two-baryon bound states 
without any additional parameters. 

There are two major factors that affect the binding behavior of the two- 
baryon system. One is the symmetry property of the system, which is charac- 
terized by the matrix-element (.A36 ), and the other is the interactions between 
quarks including both the direct and exchange components, especially the in- 
teractions caused by chiral-quark field couplings. 

We classify two-baryon bound states (even almost bound states) into three 
types. The first type state is the deuteron-like state. The symmetry structure 
of the system in this type shows that the quark exchange effect is not so 
important, namely, there is no serious Pauli blocking effect. If the chiral field 
can provide enough attraction between interacting baryons, the bound state 
can be formed. This kind of bound state might also be obtained by using the 
models on baryon level. For instance, in terms of Nijmegen model, one also 
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predicted some bound states in the 5 = — 2, — 3 and —4 systems E7\ . 



The AA-hke state is the second type state. The system of this type is 
mostly symmetric in the orbit space. In these states, the interaction, espe- 
cially the scalar meson induced interaction, dominates the binding behavior. 
Together with the symmetry behavior of the system, the relatively deeply- 
bound state can be formed. However, because both baryons in the system 
have strong decay modes, only when the predicted binding energy is lower 
than the threshold of the strong decay channel, say the NNrnr channel for the 
A A system and the ANtttt channel for the E* A system, the width of the bound 
state could be narrow. In our calculation, the predicted binding energies of 
AA and S*A are not large enough, so that the widths of these bound states 
should be rather broad. Although this kind of bound state might not easily 
be detected in experiments, it may worth to search in the future. 

The third type state is the r2f2-like state, which is the most interesting 
state in our study. Same as those in the second type, the system in this type is 
mostly symmetric in the orbit space. Meanwhile, the strange chiral fields can 
offer rather strong attraction. As a result, these states are deeply bound states. 
If we believe that the chiral SU{3) quark model is one of the most suitable 
models in describing the system with high strangeness, the predicted binding 
energy of 'E*Q is possibly below the threshold of the strong decay channel Sfivr. 
Thus, both QQ and S*r2 can merely have weak decays, and consequently are 
deeply bound states with narrow widths. 

It should be specially emphasized that the states of the second and third 
types possess the six-quark structure, and their inter-baryon distances are 
relatively short. These characteristics cannot be provided by the model on the 
baryon level. Thus, they are new dibaryon systems. The existences of this kind 
of dibaryons would be an important place to reveal QCD phenomenology. 

Appendix. 

In the framework of RGM ||2^, the upper bound is given by the expectation 
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value of the Hamiltonian 

<^) = iW- "" 

where H = T + V is the Hamiltonian, with T and V being the kinetic energy 
operator and the potential operator, respectively, \E' represents the trial wave 
function of the system and (H) denotes the upper bound of the system. The 
kinetic operator of a six-quark system can be written as: 

6 

T = Y.T,-TcM, (25) 

i=l 

where Tj and Tcm denote the kinetic energy operators of the i-th quark and 
of the center of mass motion (CM), respectively. Substituting the trial wave 
function (Eq.(21)) and the antisymmetrizer A' (Eq.(22)) into projection equa- 
tion 

{[4>A^BXCm{Rcm)]st\H - E\A'[^A^BXrel{R)XCM{RcM)]sT) = 0, (26) 

one obtains RGM equation 



{ - ^Vl+Vrt{R)-Erel}Xrel{R) (27) 

+ JiK'^iR, R') + K^'i^R, R') - EtotN'^^'^R, R')]Xrei{R') dR' = 0, 



where V^J{R) represents the relative potential between two clusters, Etot and 
Erei denote the total energy and the relative energy between two clusters, re- 
spectively, and K^{R,R'), K^{R,R') and N^''^^{R,R') describe the kinetic 
energy exchange kernel, the potential energy exchange kernel and the normal- 
ization exchange kernel, respectively. Expanding Xrei{R) by well-defined basis 
functions 

n 
Xrel{R') = T.^^^rel{R',S,), (28) 

and left- multiplying (preiiR, Si) to RGM equation and integrating over R and 
R', one obtains the secular equation of the bound state problem. 



J2[H^, - Etot N,,] c, (29) 

n 

J2[H:^-EreiN,,]c, = 0, (^ = l,---,n) 
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with Hij and A^^ being the Hamiltonian and normahzation kernels, respec- 
tively, and 

Etot = Ein + Ereh (30) 

Etot, Ein and Erei being the total, inner cluster and relative energies, respec- 
tively. Apparently, H' can be written as: 

H-j = Hij — EinNij (31) 

= [T,^ - ElN,,] + [V,, - ElN,,]. 

In this equation, Tij and Vij is the kinetic energy and potential energy ker- 
nels, respectively, and the superscripts T and V denote the kinetic energy and 
potential energy parts, respectively. Then, 

Hij' = Tij - Ei^Nij (32) 



/rpdir\ I frpexch\ 
y-'-rel Jij ~r {■'■ rel Jij 



with 



and 



(T*[)., = {<l>reliR", S,)\i-^VJiMR" - R'MreliR', S,)) (33) 



{T::i''% = {<Prel{R", Si)\K^{R", R'MreliR', S,)). (34) 



To compute Hfj' explicitly, we employ the Generating Coordinate Method 
(GCM) technique, which is equivalent to ROM. In GCM, we re-write 

n 

■^ = J2ci^i (35) 



i=l 



with 



and 



V'i = A'[(t)A<t>B<PreliR'^Si)XCMiRcM)XSTXc] (36) 

3 C. 6 c. 

= A'lYlMrn- ^)Yl(l)os{rk + ^)xsTXc] 



^os(r„-y) = Ar„exp- — (r;--)^ (37) 
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To ensure the total CM motion can explicitly be separated, we keep oj un- 
changed, namely 

u = (mX)'' = {rn,hl)-\ (38) 

If one temporarily ignores the flavor symmetry breaking, namely b = bu = bs 
and m = rriu = mg, for instance in the (QQ) or A A or A^A^ system, the 
following discussion would be flavor independent and it is easy to obtain the 
kinetic energy induced S-wave adiabatic effective interaction 



El = ^-^ (39) 

1 - 'd{P^t)Cosh{x)/Cosh{?,x) , 





3 


2Kox[ ■ 


-1 



Ctgh{?,x) 



27 {P^t)Sinh{x)/Sinh{3x) 



with 



n2 



and 



Ko = ^hu. (41) 

According to the definition of the binding energy 

Eb = -{MDB-MA-MB), (42) 

where M^^b, Ma and Mb denote the masses of the dibaryon, baryon A and 
baryon B, respectively, the contribution of the kinetic energy to the binding 
energy of the dibaryon in the adiabatic approximation should be 

n 

El = -T.cMiS.) (43) 

with 

EliS.) = El + ^, (44) 

where Cj is the wavefunction. We plot E'[{Si) versus Si in Fig. 4 (for conve- 
nience, we drop the subscript i in Si in Fig. 4). In this figure, the solid and 
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dashed curves denote the results with (-Ps^g ) = —1/9 and (-P36 ) = —1/81, 
respectively. Apparently, when (-P36 ) = —1/81, E'^{Si) has less attractive 
character than that in the (-Pa^e ) = —1/9 case. And when (P36 ) > 0, it 
would be even less attractive. 

Apparently, E^ is always a negative number (or repulsive for binding). 
This indicates that, when off-diagonal matrix elements contribution can be 
neglected under adiabatic approximation, E^ always has a repulsive effect. 

However, when calculating E^, one should also take the off-diagonal ma- 
trix elements into account. Further calculation shows that the contribution 
of the off-diagonal matrix elements differs substantially in the system with 
(-^36 ) — —1/9 to that in the system with (^35) ~ 0. In the former case, 
it is so large that E^ can be changed to a positive value (or attractive for 
binding). For instance, in the ^7^2(0,0) dibaryon, the contribution of the diago- 
nal matrix elements is —145.6 MeV while the contribution of the off-diagonal 
matrix elements is 154.2 MeV, and consequently E'[ = 8.6 MeV. On the other 
hand, in the latter case, this contribution is not large enough to change the 
sign of E^. For instance, in the deuteron case, the contribution of the diago- 
nal matrix elements is —57.5 MeV while the contribution of the off-diagonal 
matrix elements is 43.0 MeV, and consequently E'[ = —lA.bMeV. The effect 
of the off-diagonal matrix elements can also simply be seen by diagonalizing 
the —{Ejj + ■^), called Ej^, matrix. In the (P^q) = —1/9 case, some of 
diagonalized matrix elements Ef^ can be positive. Because 

n 

El = Y.c]El, (45) 

i 

where Ej^ is the diagonalized matrix element, and q is the wave function 
after the unitary transformation that diagonalizes the —{Ej, + ■^) matrix, the 
resultant E]^ becomes positive. While for the system with {P^q) ~ 0, the 
diagonalized matrix elements are always negative, consequently, E]^ always 
keeps a negative value. 

The physics picture of the positive contribution (ot attractive contribution) 
of the kinetic energy operator to the total binding energy may be understood as 
the following. In the dibaryon with {P^q) = —1/9, such as ^2^2(0,0)) six quarks 
interfere with each other very strongly due to the large quark exchange effect. 
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so that they cannot move as freely as they do in two independent baryons. 
So the expectation value of kinetic energy operator in such system (which is 
definitely positive) would be smaller than the kinetic energy in two independent 
baryons. While in system with (-Ps^e ) ~ 0, the quark exchange effect is so 
weaker that the expectation value of kinetic energy operator is always larger 
than that the kinetic energy in two baryons. 
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Figure captions 

Fig.l. A-P scattering. 

Fig.2. Energy level sketches for a. [2^7-2*^7] (1,1/2), b. [SS](o,i), c.[A^fi] (2,1/2), 

d.[Af]] (3,3/2), e.[S*A](o,5/2) and [S* A] (3,1/2), f.[AA](o,3) and [A A] (3,0), g. [f^l^] (0,0) , 

h. [S*(]] (0,1/2) • 

Fig. 3. Binding energy Eh with respect to b^ and rris for systems with {A^'^'^) ~ 

2. The solid curves are for ms=470 MeV, the dashed curves are for ms=515 

MeV. 

Fig. 4. In the S-wave case, the contribution to binding energy provided by the 

kinetic-energy related terms with different S in the adiabatic approximation. 
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